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Abstract. We investigate the cosmic reionization history by comparing semi-analytical 
models of the Lyor forest with observations of high-z quasars and gamma ray bursts absorp- 
tion spectra. In order to constrain the reionization epoch z rei , we consider two physically 
motivated scenarios in which reionization ends either early (ERM, z re j > 7) or late (LRM, 
Zrei ~ 6)- We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.7 < z em < 6.4 
and in the spectrum of the GRB 050904 at z = 6.3, studying the wide dark portions (gaps) in 
the observed absorption spectra. By comparing the statistics of these spectral features with 
our models, we conclude that current observational data do not require any sudden change 
in the ionization state of the IGM at z ~ 6, favouring indeed a highly ionized Universe at 
these epochs, as predicted by the ERM. Moreover, we test the predictions of this model 
through Lya- emitters observations, finding that the ERM provide a good fit to the evolution 
of the luminosity function of Lya emitting galaxies in the redshift range z = 5.7 - 6.5. The 
overall result points towards an extended reionization process which starts at z > 11 and 
completes at z re j > 7, in agreement with the recent WMAP5 data. 
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1. Introduction 

In the last few years a possible tension 
has been identified between WMAP5 data 
dDunklev et al. Il2009h and SDSS o bservations 
of qu asar (QSO) absorption spectra dFan et al. I 
2006), the former being consistent with an 
epoch of reionization z re ; ss 11, the latter 
suggesting z re ; « 6. Long Gamma Ray Bursts 
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(GRB) may constitute a complementary 
way to study the rei onization process, pos- 
sibly probing z > 6 ( Taglia ferri et al. I l2005t 
iGreiner et all 120091: ISalvaterra et al. I £0091) 
Moreover, an increasing number of Lyman 
Alpha Emitters (LAE) are routinely found at 
z > 6 (IStark et al. 11200% . 
The ultraviolet radiation emitted by a 
QSO/GRB can suffer resonant Lya scattering 
as it propagates through the intergalactic 
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Fig. 1. Evolution of the neutral hydrogen fraction. 
Thick red solid (blue dotted) lines represent average 
results over 100 LOS for the ERM (LRM), while 
the thin lines denote the upper and lower neutral hy- 
drogen fraction extremes in each redshi ft interval. 
Solid circles represent x H i estimates by iFan et al. I 
(2006); empty squares denote the results obtained 
in this work. 



neutral hydrogen. In this process, photons 
are removed from the line of sight (LOS) 
resulting in an attenuation of the source flux, 
the so-called Gunn-Peterson (GP) effect. For 
these reason QSO/GRB absorption spectra 
are recognized as very powerful tools for 
measuring the neutral hydrogen fraction xhi 
in the InterGalactic Medium (IGM), and 
hence for determining z le i- Moreover, the 
observed Lya Luminosity Function (LF) of 
LAEs can be used to infer the ionization state 
of the IGM at redshifts close to those of the 
star-forming galaxies and hence to reconstruct 
the reionization history. In this work, we 
compare the predictions of theoretical models 
of reionization with QSOs, GRBs, and LAEs 
observations. 

2. Reionization models 



refined by IChoudhurv et al. I (120081) . This 
model, hereafter called CF05, allows to build 
a wide range of reionization scenarios which 
match many observational constraints (ranging 
from WMAP5 to SDSS data) and differ for 
z le i and for the IGM properties at z > 6. 
In order to constrain reionization we have 
adopted the CF05 model to set up two different 
reionization scenarios, namely: (i) an early 
reionization model (ERM), favored by WMAP 
data, characterized by a highly ionized IGM 
at z ~ 6, and (ii) a late reionization model 
(LRM), in which z le i ~ 6, as suggested by 
the SDSS results. The jchi predicted by the 
two models is shown in Fig.l. Moreover, we 
have developed a semi-analytical model to 
simulate QSO absorpt ion spectra starting fro m 
the CF05 predictions dGallerani et al. 112006b . 

3. RESULTS 

3.1 . QSOs absorption spectra 

We use obse rvational data inc luding 17 QSOs 
obtained by IFan et al. I (120061) . We divide the 
observed spectra into two redshift-selected 
sub-samples: the "Low-Redshift" (LR) sample 
(5.7 < z em < 6), and the "High-Redshift" (HR) 
one (6 < z em < 6.4). By comparing the the 
Largest Gaj0 Width (LGW) distribution ob- 
tained from the observed spectra with simula- 
tions (Fig. 2) we find that both the ERM and 
the LRM provide a good fit to observational 
data. We exploit this agreement to derive an es- 
timate of jchi. We find log 10 XHi = -4.4*9gi at 
z ~ 5.3, and log 1() jc ffl = -4.2^;jj 4 at z ~ 5.6. 
Although the predicted LGW distributions are 
quite similar for the two models considered, in 
the HR case we find that a neutral hydrogen 
fraction at z ~ 6 higher than that one predicted 
by the LRM would imply an even worst agree- 
ment with observations, since a more abundant 
HI would produce a lower (higher) fraction of 
LOS characterized by the largest gap smaller 



We present a semi-analytical model 

of cosmic reionizatio n, developed 

by IChoudhurv & Ferrara 
Choudhurv & Ferrara 



(120051) and 
(120061) and further 



1 Gaps are defined as contiguous regions of the 
spectrum characterized by a transmitted flux lower 
than a flux treshold F t h = 0. 1 . 

2 The LGW distribution quantifies the fraction of 
LOS which are characterized by the largest gap of a 
given width. 
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(higher) than 40 A with respect to observa- 
tions. Th us, this study suggests xhi < 0.36 at 
z = 6.32 dGallerani et al. Il2008al) . 

3.2. GRBs absorption spectra 

We have analyzed the optical afterglow spec- 
trum of the GRB 0509004 at z = 6.29 by mea- 
suring the LGWs in its Lya forest as a func- 
tion of the flux threshold f\h used to define 
gaps. Starting from a large sample of synthetic 
spectra, we have computed the probability to 
find gaps of these widths in absorption spectra 
obtained through the ERM and LRM (Fig. 3). 
We find that the ERM is two times more pre- 
dictable than the LRM, thus confirming the re- 
sults found in the case of QSOs, i.e. that a 
highly ionized Unive rse at z ~ 6 is favored b y 
current observations (Gallera ni et al. 1 12008b). 

3.3. LAEs luminosity function 

We have built up a semi-analytic model of 
LAEs to compute the evolution of the Lya 
LF of star forming galaxies at redshifts ap- 
proaching Zrei- We have fixed x H i to the pre- 



dictions of the ERM and LRM at z = 5.7 and 
z - 6.6 leaving only the star formation effi- 
ciency and the effective escape fraction of pho- 
tons as free parameters. The results are shown 
in Fig. 4. The ERM provides a good fit to obser- 
vational data, implying that no sudden change 
in the ionization state of the IGM at z ~ 6 is 
requ ired to explain the evolution of the Lya 
LF (iDaval et al. 1 120081) The ER M has been 
also used bv lDaval et alH d2009al) to predict the 
LAEs L F up to z — 7.6 and by IDaval et al. I 
d2009bl) to study the dust content of LAEs. 



4. Conclusions 

We have presented a semi-analytical model 
which allows to simulate QSOs, GRBs, and 
LAEs spectra. Synthetic spectra have been an- 
alyzed statistically and the theoretical predic- 
tions compared with observations. We have 
shown that current data do not require any 
sudden change in the IGM ionization level at 
z ~ 6, favouring a highly ionized IGM at these 
epochs. The overall result points towards an 
extended reionization process which starts at 
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Fig. 3. Isocontours of the probability that the afterglow spectrum associated with a GRB at redshift zgrb = 
6.3, contains a largest gap of size in the range [W max , W max + dW], with dW = 20 A, for a flux threshold F th . 
The left (right) panel shows the results for the ERM (LRM). The isocontours correspond to probability of 
5%, 10%, and 40%. The black points indicates the position in the (W max , F t ° h bs ) plane of GRB 050904. The 
point with arrow means that the gap size should be considered as an upper limit, since the corresponding 
dark region could be affected by the presence of a DLA. 




Fig. 4. LAE LF for the ERM at z = 5.7 (dot- 
dashed), z = 6.5 6 (solid). Points represen t the data at 
z = 5.7 (circles. lShimasaku et aTl d2006h ). z = 6.56 
(down ward/upward triangles re present upper/lower 
limits; iKashikawa et al. I O006h ). 
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z > 1 1 and completes at z re ; > 7, in agreement 
with the recent WMAP5 data. 



